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THE MOTIONS OF ROLLING SYMMETRICAL MISSILES 
REFERRED TO A BODY-AXIS SYSTEM 
By Robert L. Nelson 


SUMMARY 


The linearized equations of motion have been derived for a rolling 
missile having slight aerodynamic asymmetries. Time histories of rolling- 
missile motions referred to a body-axis system, have been prepared to show 
the types of missile motions that can be encountered. The motions 
resulting from a trim change and a pulse -rocket disturbance are shewn to 
be determined mainly by the ratio of rolling velocity to pitching frequency. 

Finally, the derived equations are used in establishing a technique for 
the reduction of rolling -missile oscillation data. It is shown that the 
aerodynamic derivatives can be obtained from flight data if four accelera- 
tions are measured. The method is applied to the results obtained from a 
flight test of a missile configuration. 


INTRODUCTION 


The problem of the dynamic stability of missiles has been attacked by 
two separate basic treatments. First, in the case of roll-stabilized mis- 
siles the problem has been attacked by means of the classical airplane- 
stability theory developed by Lanchester (ref. 1), with no need for modifi- 
cations. The development of the theory for rolling symmetrical missiles 
has followed from the basic, ballistic theory of reference 2. Nicolaides 
(ref. 3) and Charters (ref. 4) have more recently expreshed the dynamic- 
stability equations for rolling symmetrical missiles in terminology more 
familiar to the aerodynamicist . However, the impetus for these works was 
provided by studies of missile motions in aerodynamic test ranges where the 
model position and angular orientations were measured at a series of sta- 
tions. As a result the equations derived were referred to space axes. In 
the case of free -flight rolling missiles equipped with internal instrumen- 
tation, the equations of motion must be referred to a body-axis system. 
Phillips in reference 5 presented a simplified analysis of the motion of 
rolling airplanes (and in a limiting case of symmetrical missiles) based on 
a body-axis system but included only the criteria for stability. Bolz in. 
reference 6 derived more completely the equations' of motion for a rolling 
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symmetrical missile referred to a "body-axis system "but again discussed 
only the criteria for stability. 

In this paper the equations of motion referenced to the hody-axis 
system are again derived, and assumptions similar to those of Kicolaides 
(ref. 3) are used. Same of the possible missile motions are shown, 
together with the motions to be expected from certain forcing functions. 
Finally, the derived equations of motion are used to establish a techni- 
que for the reduction of oscillation data to obtain aerodynamic deriva- 
tives. The method is applied to the experimental results obtained from 
a rolling symmetrical missile. 
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SYMBOLS 

resultant acceleration 

complex roots of differential equation for £ 

accelerations parallel to X, Y, and Z axes, 
respectively 

constants of differential equation for £ 
body frontal area, sq ft 

My 


pitching-moment coefficient. 


qApd 


Ss - I s ’ “■= 
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"N 


normal -force coefficient (-Cg) 
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c n 


yawing-moment coefficient, 


Mz 

QAjid 




' n P 

Bp 


■ 


Bijr * 

1 

_ SCn 


sec 


^ a " 

, _ ^Cn 


\ 
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C F magnitude of resultant vector 

c R,o = % + ^2 


Cy 

^,0 


C = 

Y P Bp 


C Z 

C Z,o 

n = ^Z 

^ Bo. 


fy 

lateral- force coefficient. — — 

q.A F 

lateral-force coefficient due to asymmetry 


vertical-force coefficient, — — 

qA F 

vertical-force coefficient due to asymmetry 


d body diameter, ft 

F X> F Y> f Z aerodynamic forces parallel to X, Y, and Z axes, 
respectively 

I moment of inertia about Y- or Z-axis, slug-ft^ 


I* = 


qA F d 


sec c 
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k 


moment of inertia about X-axis, slug-ft^ 


K - R X-~Jg. 
R 1 + ®2 

k 

l 

My, Mg 


radius of gyration, ^i/m, ft 

pulse -rocket moment arm from center of gravity , ft 
aerodynamic moments at out X, Y, and Z axes, respectively 


m 


mass, slug 


m' = 


mu 

qAj.. 


sec 


P 

P 

q. 

R x , Rg 


pulse-rocket total impul.se, lb-sec 
ro llin g velocity, radians/sec 
dynamic pressure, lb/sq. ft 
initial magnitude of rotating vectors 


s 

t 

u 


V 

X, Y, Z 
x o 

*1, X 2. 
?i> y 2 
* 1 , z 2 


slope 

time 

component of free-stream velocity in X-direction, 
ft/sec (u V) 

free-stream velocity, ft/sec 
coordinate axes 

longitudinal center of pressure of asymmetric force, ft 
longitudinal displacement of accelerometers, ft 
lateral displacement of accelerometers, ft 
vertical displacement of accelerometers, ft 
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a 

P 

£ = 0 + ia 

t 

e 

A 


angle of attack, radians 
angle of sideslip, radians 

pitching velocity, radians /sec 

a n gula r orientation of resultant vector, radians 


Aq = = vg, radians 


A 0 nonrolling damping constant, l/sec 

damping constant due to roll, l/sec 
M- relative -density factor, _Ji2L_ 

pApd 

v l* v 2 initial an g ul a r orientation of rotating vectors, radians 


P air density, slugs/cji ft 

• 

0 rolling velocity, radians/sec 

t yawing velocity, radians/sec 

% "basic oscillation frequency, radians/sec 

Ao) component of total pitch frequency resulting directly 

from roll, radians/sec 

ft = 0 + iijf 

Dots over symbols indicate time derivatives. 


STATEMENT OF THE PROBLEM 


The motion of a symmetrical missile configuration referred to a "body- 
axis system is co mp l i cated "by rolling motion even for the case where the 
rolling velocity is a fraction of the missile pitching frequency. Unlike 
the motion of a slowly rolling airplane -like configuration which .tends to 
continue in a direction normal to the wing plane, the motion of a symmetri- 
cal missile configuration disturbed in one plane 'in space tends to remain 
in one pla n e in space. (The restoring force for the airplane is normal to 
the plane of the wings, whereas the restoring force for the missile is 
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mainly independent of roll attitude.) Thus for a missile having acceler- 
ometers mounted at the center of gravity the first consequence of the 
rolling motion is the appearance of a roll effect on the normal and 
transverse accelerometer traces. An . example of this effect for the normal- 
accelerometer trace is shown in the following sketch: 



This roll effect is determined mainly by the ratio of the rolling 
velocity to the oscillation frequency, or the ratio of pd/2V to a^d/2V. 

At very low values of pd/2V corresponding to very small fin misaline- 
ments, this ratio can become large if the missile ascends to a high alti- 
tude, inasmuch as • pd/2V is independent of altitude while u^dyfev 

decreases with increasing altitude. 

If the missile motion were in one plane in space and the roll rate 
small, the motion could be analyzed simply by working with the resultant 
of normal and transverse acceleration, and applying the us ual data- 
reduction techniques for oscillating nonrolling models (ref. 7 ). However, 
since the missile may not remain in one plane In space and gyroscopic 
cro s 8 -coupling can have a large effect on the model motion, an analysis 
of the motions of rolling missiles is necessary in order to interpret 
and reduce the flight test data. 


DEVELOPMENT OF EQUATIONS OF MOTION 


The linearized equations of motion for a symmetrical missile (90° 
rotational symmetry) are referred to the body-axis system of the following 
sketch: 
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On the assumption that 
that large rolling velocity 
reduce to 


ang^ar deviations from trim are small hut 
must he considered, the equations of motion 


ni 8 ^ — 

* 


mu(p + f - $a) = Ef y z/= coviApoueui of -free-Jf/ean, 

Velocity ill X- di\ScWw^ 


mu(& - 9 + $p) = ZF Z 

ij = 2Mx 

19 - (i - l X )jZ|£ = 2^ 


It + (i - Ix)09 = ZM Z 
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If the forward velocity is assumed to "be constant (ay = 0) and the 
rolling velocity also to he constant (0 = 0), Idle equations reduce to the 
following four: 


mu(p + t - 0a) = ZFy 

(1) 

mu(d -9 + 00) = 

(2) 

10* - (I - l x )0t = 2% 

(3) 

It + (I - I x )00 = m z 

00 


In order to reduce the number of equations to two, equations (2) and (4) 
are multiplied hy i = yj- 1 and added to equations (l) and (3), respec- 
tively, to obtain 


rau[(p + id.) - i(0 + it) + i$(0 + iaj] = £(Fy + iF z ) 


1(0 + it) + i0(l - I X )(0 + it) = 2 (My + 1M Z ) 


Let £ = 0 + ia, a = 0 + it, 

t + iPS 


. / I x \ 2 (My + iMg;) 

n + ip(i - -f)a I , (6) 

Allowing small asymmetries (resulting from a fixed control deflection for 
example) which do not violate the assumption of symmetry and assuming that 
the aerodynamic forces due to t, 0, and 0, and forces due to gravity 
are negligible in comparison to the forces due to angle of attack and 


and 0 = p so that these equations become: 
. m = (5) 


mu 
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sideslip permits the force equations to he written as: 

* 


ZFy = Cy^qApP + Cy qAp 


ZF Z = Cg^qApa + C z ^qAp 


Then, multiplying the second equation hy i and adding to the first gives 

S(Fy + iF z ) = qAp^Cypp + iC Z a °^ + (pY 0 + ic Z Q ) 

From the assumption of symmetry, 

% = °Zd = tt 

The derivative rather than Cy^ or 1 b used hereinafter, since 

the force due to a positive angle of attack is usually considered positive. 
Then, 


Z(Fy + iF z ) = -C|j a qAp(p + ia) + qAp ^Cy Q + iC Z( ^ 


Let m' = then, inasmuch as £ = p + ia, 

qAp 


Z(Fy + iF z ) = _ + C Y 0 + iC Z 0 

mu m' m' 


L 
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If no moments from asymmetries . are e allowed to result from camber 
effects and if the moments due to 0, tJTj and 0 are assumed to mak e a 
significant contribution to the total moment, the moment equations can 
he written as: 


= C^qApdo, + C^qA-pdfi + C^qApd^ _ C^qA*. 

£Mz = C np qApdp + C^qApCbfr + qApd&t + Cy o X 0 qA F 

Again, multiplying the second of these equations by i and adding it to 
the first gives 

E(% + 1M Z ) = qApdHCi^a + + (Cj^0 + + 

- 10 ^) " K " 1Ci o)£ 


With the assumption of 90° rotational symmetry. 


% " "Sc, 

S = S 

c % = % 


♦ 

J 
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Then 


S (My + tM^) - -igApd + 5^(0 + ^ a ) + 

10^(0 + it) - (Cj 0 + ic Zo )^J 

Also, with £ = 0 + ia, ft = 8 + it, $ =■ p, and I* = - J --, 

QApa 

= ' r ( c “a + ipCn &^ + iCm e a " (% + 1 c z 0 )t 


Equations (5), (6), ( 7 ), and (8) combine to give two linear differential 
equations in terms of the two unkno wns £ and ft. The following equa- 
tion in' £ results: 


£ + A£ - B£ = C 


where 




nr I 1 


<■-*) 


= i C 3sk p + p 2 A . a\ + . lp & _ _ 5c 5Ss\ 

. I* \ 1/ I' m’ \m’ I* I m'/ 

- - " C 

0 „ . K . !3l . ipA . is) (% Iggol (: 

— 
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The solution of the differential equation is: 


{ = %* 1Vl * at - 


-A + yllT + 4B 


-A - VA 2 + hB 




iv l lv 2 

and R-j^e and Rge are initial conditions as shown in the adjacent 
sketch. 



Both a and h are complex numbers , with the complex radical 
/a 2 + 4B, and can be broken down into a more usable form to 
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a = Aq + AX + i(flct) - Am) 
"b = Aq - AA - i(flCb + Aa>) 


where the new concepts Aq, AA, and Ao are defined by 



AO) = p^L - e. I (l6) 


From the new equations for a and b, the following expression can he 
obtained: 



f = . B e iv l e [^D+AM-i (%-Ao>|'t + Roe iv 2e [Ao-AA-iCo^+Am)]! 

•* •’t rim 1 *2 

or 

^ = ^ + Rie (Ao + AA)t e i[(a^nAaj)ti-vJ- + 

(17) 
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Inasmuch as a complex number represents a vector which can he expressed 
in the form z = a + ib = Re i0 , where R represents the t magnitude of the 

ma 


vector and 8 gives the angular orientation. 


R 2 e 


(^o^A)t 


R^e' 


represent the damping with time of R-j_ and R 2 , respectively 


and (ob - Aoj)t + V]_ and (cdq + Aco)t - v 2 give the location of the vec- 
tors at any time. These relationships axe shown in the following sketch: 



Thus, the resulting motion is given "by two clamping vectors revolving in 
different directions and at different rates about the fixed trim. 

The time histories of a and 0 can he obtained from equation (17) 
by using the relations e 10 = cos 0 + i sin 0 and t; = 0 + ia; then 

(P " Ptrim) + 1 ( a - a trim) 


cos [(cut, - Aa>)t + vj + 

_ _ 

Rge cos Kab + - v 2 j + 

sin|(ab - £u>)t + vj| - 

0\ o -£\)t 

R2 e sin (at, + Ao)t - v 2 l 


( 18 ) 
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and 


(V^x)t 

P - P-trim = R l e cos [fao - + v L J + 

(Ao*AA)t p. 

R^ cos Ho^ + An)t - V2J 


|— . -1 

a " O'trim " E l e sin [fab - + vjj - 

(Aq"'^)'^ I -. . -j 

H2e sin (cdq + Aa>)t - v 2 


( 19 ) 


Thus the oscillation about trim of either a or p is made up of two 
sinusoidal components which can have different initial magnitudes 
(Rl £ R2) and which can damp at different rates (AA £ 0) . When AA 
approaches the value of Aq, one of the component oscillations damps 
rapidly while the other component damps slowly. 

Rather than use equations (l8) and (19) in the discussion which 
follows, use will he made of the resultant angle of attack from trim and 
the angular orientation of this resultant. From equations (18) and ( 19 ) 


(5 - ♦ (• - ♦ 


2R l R 2 e2 ^° t cos(2obt + v 1 - v 2 ^ (20) 

a - = _ &p)t + yj - sln[(ab + am)t - vgj 

^ ^trim ^ cos J(a^ + am)t - vgj 


One interesting point should he made concerning equation ( 20 ). The roll 
rate p appears as an important quantity only in the quantity AA, which 
is small compared with Aq at lew roll rates. Thus, to the first order, 
a change in the rolling velocity will not influence the resultant 
oscillation time history, hut will only influence the angular orientation 
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of the resul tant . Therefore, the assumption that the rolling velocity is 
constant may not he too strict a limitation on the validity of the >• 

analysis at low roll rates. 

T’hft sections which follow discuss the motions of rolling missiles ' 

referred to the body-axis system and a technique for reduction of data of 1 
missiles employing body-axis instrumentation. In this system forces 
measured at the center of gravity more conveniently describe the motion 
t han would, the angle of attack, which cannot be measured at the center of 
gravity and must be corrected to the center of gravity. 

From the assumption mndp earlier that aerodynamic forces due to 0, 
ijr, ^ are negligible, the equations describing the force time his- 

tories differ from derived equations only by a constant. If the forces 
due to 8, \jr, and $ are not negligible, the new equations describing 
the force time histories differ from equation ( 17 ) by the initial condi- 
tions R-l, Rjjj v-l, and v 2 - 

In the equations and figures to follow, use is made of the symbols 
as defined in the following table, for either angle -of -attack and side- 
slip instrumentation or for force -measuring instrumentation. The coeffi- 
cient C z , which is positive downward, has been replaced by -Cjj. 

• t 1 


Symbol 

a and (3 system 

Cjj and Cy syatem 

°R 

\/( P " Ptrim) + ( a “ “trim) 

\jipY " ^trim) 2 + (°N - C N,trim) 2 

pR,trim 

/ 2 2 
y “trim + Ptrim 


A 

tan-1 “ “ afcrim 

P ~ Ptrim 

« _i °N “ c N,trim 
it - tan x 

“ °Y,trim 

Afcrim 

tan" 1 atrim 
Ptrim 

_1 ^.trim ” ^,0 
it - tan x —*■ — 2 — - 

c Y,trim “ °Y,o 
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DISCUSSION OF THE DERIVED MOTION 
/ 

In the discussion to follow, motion is intended to mean the missile 
motion as reflected hy the resultant angle-of -attack or resultant force 
time histories. 


Rolling Trim 

By use of the equation tw| = -C/B, the rolling trim referenced to 
the nonrolling trim can "be expressed as : 
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The quantities N and Ajj are independent of the stability roots \j, ( 

£D\, and Zko. ‘In addition, N and Ajj are' independent of velocity, 

and independent of Mach number (if pd/2V, C^, and Xo/d are constant) , 

and N is influenced slightly by altitude changes (through the 
CmqfV- term) . 

Same examples of rolling-trim magnitude and orientation are given in 
figures 1 and 2 for a particular missile configuration. The assumed 
aerodynamic and mass characteristics listed in the figures are representa- 
tive of a four-fin missile of high fineness ratio. Tha resonant condition 
at Ao/ab = tl is shewn in figure 1. An indication of the magnitude of 

the rolling trim at resonance for low values of damping is also shown. 

Figure 2 shows that the trim rotates through l80° as the resonance range 
is passed. 

In order to give a better idea of the importance of the resonance range ■ 
and the associated problem of avoiding resonance, the equations for 'h o /(n 0 
and Aoi/cub are first written as follows: 



(26) 


• « 


( 27 ) 
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Note that both equations ( 26 ) and ( 27 ) are functions of only the mass 
characteristics of the missile, the aerodynamic characteristics, and the 
altitude (through the relative density factor) and are independent of 
velocity, except indirectly through the effects of Mach number and 
Reynolds number on the aerodynamic derivatives. 

Figure 3 shows the variation of and &a/ act, wit h altitude 

for a high-speed missile . The asymmetries causing the roll were acci- 
dental and probably on the order of what can be expected with quality 
construction of a finned missile. Ohe curves show resonance reached at 
an altitude of 93 >000 feet and at 'K q ^(Tq = -0.022. Passing through 

resonance at such a high altitude may be objectionable; for, even though 
the resultant force may be less than at a lower altitude, the greatly 
increased trim angle of attack may be such sis to put the missile in an 
unstable angle -of -attack range, which could lead to destruction of the 
missile or at least to a significant change in flight path. In addition, 
if the flight path is such that resonance occurs near zenith, then the 
resonance range will be passed through slowly and large changes in flight 
path can occur. It may be that, for certain applications, the fins should 
be deflected, so that resonance occurs at a predetermined low altitude. 
Although the motion about trim is not greatly affected by changes in 
rolling velocity at low roll rates, this result does not hold true for 
the trim if changes in roll i ng velocity take place near roll resonance. 


Motion About Trim 

Setting zero time at an, oscillation peak (v^ = v 2 ) and replacing in 
equations (20) and (21) the initial conditions R^, Rg, v^, and v 2 by 
the following: 

c R,o “ R 1 + **2 


Aq = v L = v 2 

V ' K . 

' R]_ + Rg 

yield the following simplified equations: 

= + K ) ?e2AAt +' C 1 - K) 2 e -2AXt + 2(1 - k2)cos 20^ ( 28 ) 

A 

% 

\ ( 


l 
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A = Aq 


- £ajt + tan -1 


(1 + K)e^^ - (1 - K)e^ 
(1 + Kje^ + (1 - K)e^ A?vt 


tan afct 


(29) 


A "better insight into the physical significance of the quantity K is 
obtained from the equation for dA/dt: 


^ _ cob [(I + K^e 2 ^ 1 - (1 - K^e" 2 ^] + 2(1 - K?)A\ sin 2ofet _ ^ 
dt (1 + K^e 2 ^ + (1 - K) 2 e _2AXt + 2(1 - K^cos 2ffie>t 

(30) 


At t = 0 


Q) = Kob - (31) 

dt /t=0 


Thp value of K is thus related to dA/dt at t = 0. New the motion 
about -trim referred to an axis system in space, which translates, pitches, 
and yaws with the body-axis system but does not roll, is related to the 
motion with respect to a body-axis system by the following equations : 



■^space -^body axes + 


since the body-axis system rotates at the roll rate p. Thus, the motion 
in space is determined by the initial conditions Ao and dA/dt at 
t = 0. Also, 



K = 1 ^space^ 
“b\ at J^q 


For values of 
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Also, for tlie condition of zero roll rate (p = 0), 



Thus, from this, the motion from trim in space is nearly independent of 
the roll rate at low roll rates and. low values of Iy/I. Since the 

motion in space is determined mainly by the value of K, the quantity K 
is called the "space -motion factor". 

Figures 4 and 5 present some typical motions about trim (Cp against 
A with time varied) given by equations (28) and (29) for variable values 
of K and Aoyh^ and for zero damping (Xq and AX = 0). The effect of 
Xq on the patterns is mainly to reduce succeeding peaks. The effect of 

AX is discussed in a subsequent section. The number of cycles to com- 
plete the pattern are also shown. 

Seme of the important characteristics of the motions shown in fig- 
ures 4 and 5 are as follows: 

(1) For Acu/ufc, = 0, or for the motion in space, the pattern is 
elliptical, varying from a line (motion in one plane) for K = 0 to a 
circle for K = ±1. 

(2) Whether K is positive or negative can be easily detected, 
for Cp circles the trim center for negative values of K, but does not 
for positive values of K. 

(3) An indication of the approach of resonance is shown by the 
tendency for the loops to close on each other and become coincident for 
a given value of K. 

(4) Whether the missile is above or below roll resonance is Imme- 
diately determined, since the motion above resonance is characterized by 
inside loops, and the motion below roll resonance by outside loops. 

(5) In order to determine the value of the rolling trim from a plot 
of Cj[ against Cy, a low value of K is desirable. 

As indicated, a plot of Cjj against Cy (or C B against 
yield immediately valuable information about the motion. 


A) will 
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Rolling influences the damped motion through the term AX. From 
equation ( 28 ) the envelopes of Cp are: 


+ K)e Mt + (1 - K)e"^^j (32) 

Cp r, 2 -* 


C R,min 

°R,o 



Kje^ 



(35) 


Typical oscillation envelopes of Cr helow divergence (AX< \ ) ) are 


shown in figure 6. The figure shows that the effect of an increase in AX 
tends to damp the oscillation of Cr more highly, except for positive 
values of K. For positive values of K the oscillation becomes less 
damped with an increase in AX/Xq, up to the time where = 0. But 


with further increases in time, the oscillation becomes highly damped. 
Whether AX/Xq is positive or negative can be obtained from the following 

equation for AX/Xq : 


AX = 


21 


C N + 
“a 




i-3 

21 


f c w a - 


«r 


Aai _ 

Cm q\ <D ° 


>fkV 

-\a.) 


S S a 


1 - 


21 


f C N a ' 


2 t 




An 


(34) 


For conventional, missile configurations 

of greater magnitude than Cjj a ; as a result 
sign as Acn/cq,. 


1 



AX/Xq 



is negative and 


usually has the same 


Since the limiting magnitude of Au/oq, is 



can be between ±1, 


AX/Xq will 



never exceed ±1, 


if 
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Cnp^ is negative. However, a positive value of Cn^ can lead to 

instability at high roll rates or high altitudes, or both. However, the 
instability resulting from a positive Cnj^ will be serious only if it 

occurs at altitudes where significant dynamic pressure exists, since the 
rate of divergence is negligible at extremely high altitudes. Thus , 
knowledge of the value of is most important for missiles which 

roll rapidly. 

Another way of visualizing the effect of AX on the motion is 
obtained by finding the variation of K from peak to peak of the 
oscillations or effectively setting the initial conditions and zero 
time at succeeding peaks. Then, from equations (30) and ( 31 ), 


K t = (1 + K)e^ - (1 - Kle^ 
(1 + + (1 - 


Since e^’ t or — *-0 as t— *-<», K-fc— >-±l as- t— Thus, AX 

causes the motion to tend toward that given by K = ±1 in figure 5 . 
Figure 6 is in agreement with this concept. Of course, at low roll rates 
this trend is completely obscured since the oscillation is, for all prac- 
tical purposes, damped out by the time the influence of AX is felt. 


RESPONSE OF ROLLING MISSILES TO FORCING FUNCTIONS 


Thus far, motions resulting from variable initial conditions have 
been discussed. Remaining to be determined are the initial conditions 
(and as a result the motions) that can be expected from given disturbances. 


Response to a Trim Change 

Assume that, at the time of the instantaneous control deflection, the 
missile was at the initial trim and the oscillation about this initial trim 
had completely damped out. The initial conditions are referred to an axis 
system at a - a trim , (3 - P-trim ® bou ' t which the oscillation takes place. 
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as indicated on the following sketch: 



£t=0 ~ 0 


vhere 


^R,trim “ y(%,trim ) 2 + (^trim^ - ^,trimfe,trim)o cos ^rim 
and 


A_ = tan ~~*- S ^ n •^trim^o ~ ^R,trim s ^ n -^trim 

(^Rjtrimjo cos Atrimjo " ^trim cos ^trim 

Substituting these initial conditions into equation (17) gives the 
following values for R-]_, Rg, and v 2 - v-j_: 
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\ ftob + Am) 2 + (\) - AX) 2 


^ R > trim 2^7^ 


R2 - &n) 2 + (\, + AA)‘ 

2(/^ 2 + AA 2 


_1 (o^ - Am) _i (% + Am) 

v 2 - v i = it - tan —r^ — - tan 


~(Ao +.• AA) 


-(Ao - AA) 


Neglecting AA gives as the maximum value of C R , from equation (20), 


c R,max “ ®*° t *(»l + ? 2 ) 
at the time t*, which is given by 


t» = V 2 H i 

2a>o 

Assume that A Q and AA sire small quantities compared with % and Am. 
For Am < m^. 


t* = 0 


C R,max - ^R,trim 


For Am > cDq, 


t* = _2L_ 
2ob 
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= ^ + ^ ryi e It V 2a5 ° 


( ^R,max “ a^j ^ c R,trim e 


The value for K is 


k „2liI2 

Rl + Rg 


Then, for An> < 


K = 

tcu 


and for Am > 

K Am/a*, 

These results are plotted in figure 7» It should not "be inferred from 
figure 7 that the value of C- Rj TnHV increases with an increase in Am/o^ 

above Am/ccfc, = 1 , since Or is referenced to ^C^trim' whi ch drops 

off rapidly with an increase in Am/me> above Am/ccb = 1 , as is shown in 
figure 1. Note that, since the space-motion factor K and Am/o^-, are of 

the same sign, only the motion patterns for positive values of K in fig- 
ures 4 and 5 are. possible. 


Response to Pulse -Rocket Disturbance 

A free-flight missile can be conveniently disturbed from its equi- 
librium conditions by filing a sttibI 1 rocket charge mounted normal to the 
longitudinal axis and forward or rearward of the center of gravity. Such 
a rocket charge is referred to as a ' 'pulse rocket." For this simplified 
analysis the rocket thrust is assumed to be constant. However the pulse- 
rocket total impulse (area unde r the curve of thrust plotted against time) 
and the pulse-rocket burning time, rather than pulse -rocket thrust, are 
used as variables . By use of methods similar to those used in the section 
"Development of Equations of Motion"", solutions have been obtained for the 
peak value of C R and the value of K at the peak, for cases where the 

peak C R occurs before rocket burnout and after rocket burnout. The 
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results are presented in figures 8 and 9 for the special case of zero 
damping (Ao = 0, AA = 0) . 

Figure 8 shows the large effect the roll rate can have on the peak 
value of C R , where C R is indicative of the resultant angle of attack. 
Figure 9 shows that the space-motion factor K rapidly approaches 1 at 
a given pulse-rocket burning time, as Acd/ouo increases in magnitude. 
Both figures show the desirability of a very low pulse-rocket burning 
time, first in order to obtain a high peak Cr occurring after pulse-, 
rocket burnout, and second in order to restrict the motion to nearly 
one plane in space. 


APPLICATION OF EQUATIONS OF MOTION TO ESTABLISHMENT OF A 
TECHNIQUE FOR REDUCTION OF OSCILLATION DATA 


The procedure used in applying the equations of motion for the purpose 
of reducing the oscillation data is presented in two parts. The first part 
deals with the determination of the stability roots, Ao, AA, o^, and 

Ao), and the second part deals with a method for obtaining the static and 
dynamic aerodynamic derivatives from the stability roots. It is shown that 
the derivatives can be found by measuring four quantities in addition to 
airspeed and atmospheric conditions. Of course, the accuracy of the deriv- 
atives obtained depends, to a large extent, on how well the assumptions 
set down in the previous analysis are satisfied. 


Determination of Stability Roots 

Consider a rolling missile employing instrumentation at the center 
of gravity for measuring the normal-force and side-force coefficients. 

If the missile is disturbed in such a manner that the space-motion fac- 
tor K is near zero and is out of the region of roll resonance so that 
the trim does not wander in magnitude and direction, plots of Cjj against 

Cy (or Cr against A) similar to those shown in figure 4 and 5 can be 
obtained. A typical experimental plot is shown in figure 10. With the 
low values of K, a fairly reliable trim center can be assumed. 

The quantity Acu/a^, is readily found by measuring the angle between 
adjacent peaks on the plot of Cr against Cy. (See following sketch.) 
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A x = Aq - 

Ag = Aq - Aujtg + 11 


Ag - A]_ = -Ao(tg - tj_) + it 


and since t 2 and are measured at adjacent oscillation peaks 


tg - % 


it 

«b 


and 


4 s = 1 - fkzii (35) 

m O 11 


Calculation of &d/o£> for all adjacent oscillation peaks gives an indica- 
tion of the validity of the assumption that the roll rate p is constant. 
Also the accuracy of the calculated Aa> is affected by the correctness of 
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the assumed trim. The following table presents the calculated values of 
Atn/cOo for the experimental plot of figure 10: 


Peaks 

Aos/flCb 

1 to 2 

0.301 

2 to 3 

.325 

3 to 4 

.282 

4 to 5 

.308 


The smal l variation of Aas/oso upholds the assumption of a nearly con- 
stant roll rate. 

With the trim center known, the time history of C R 2 can be obtained 
from the plot of Cjj against Cy* If sufficient oscillations are present, 

the damping envelopes about the oscillation peaks can be drawn with good 
accuracy. Equations ( 32 ) and ( 33 ) represent the time histories of the 
damping envelopes. If use is made of the damping envelopes, equation ( 28 ) 
can be written as: 

°B 2 - | ! c R,ma4 2 + + c R,o 2 ( 1 - K^e 27 ^ cos (56) 

Thus from this equation, C R 2 is a sinusoidal damped oscillation about a 
trim given by ||i,max) 2 + (%min) 2 | . From the period of oscillation cct> 
can be obtained. From equations ( 32 ) and ( 35 ) and the damping envelopes, 

\ log j( C R,max) ” ( c R,min)^j = ^o^ + Constant (37) 

\ log ffi^ R^min \ _ + constant (38) 

\“R,max " c R,min/ 

Thus, 7 \q and AA are obtained from the slopes of plots of 
| l0 S [Kmax) 2 - (0R,min) 2 ] - ^ 
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Figure 11 is a plot of Or 2 obtained from figure 10. Except near 

the end of the oscillation, where an error in the assumed trim would have 
a large effect on time for the oscillation peak, the periods measured at 
oscillation pe aks appear quite good. If an average period of 0.280 second 
is used. 


ab = 


it 

0.280 


11.20 


Two special corrections were also used in evaluating the damping ~h c 
First, since the assumed trim may not be correct, the error would appear 
as an oscillation in 




As a result, the periods 

measured at points where the oscillation cuts i- 

R,max ) 2 + (pR,min)j 
would be inconsistent and differ from the periods measured at the peaks . 
An indication of the trim error is given by this inconsistency of the 

period. As a result an adjusted value of ^ [( C R,max) 2 + (pR,min) 2 j *«» 

been determined for each half cycle, such that the period is consistent 
with the adjacent peak periods. These adjustments to 


| K C R,max) 2 + (^min) 2 ] 8X6 ideated by the dashed lines. 


While this 


correction is admittedly crude, it will provide a more nearly correct 
damping curve. 


The second correction is in consequence of the high speed and the 
flight path of this particular model. As a result, the air-density 
change during this plotted oscillation was significant. On the assumption 
that the change in air density does not force the model oscillation, the 
damping at any time during the oscillation (if the velocity c han ge is 
assumed negligible) can be written as 


\>,t = *o,t=0 

Then from equation (57) 

2 log [[ C R,inax) 2 " ( c R,min)^ = A o,t=0^^j 
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Thus, this correction appears as a time adjustment. Figure 12 presents 
the damping determined, including the effects of each correction. From 
the slope of the lower plot in figure 12, 


= - 0-736 

Because of the low roll rate, AX is small and cannot be obtained from 
figure 11, since (c^min) 2 0. 

Determination of Aerodynamic Derivatives 
From Stability Roots 

First, equations are written for the differential accelerations at 
two points on the missile as measured by normal, and transverse 
accelerometers : 


Aag = a2;(x2,y2,Z2) - agCx^y^zi) 


-■ - ( x 2 ~ x l) j^~ ^ - ( z 2 - Z 1 )? 


Aa y = ay(x 2 ,y 2 ,Z 2 ) - ar(»L>yi>Bi) 


= (xg - + ^ ^ ■ ( y 2 - yi )^ 2 


These equations were derived under the same assumptions that were made in 
the section "Development of Equations of Motion". Since ACjj = -mAag/qAp 

a nd ACy = mAay/qAjp, 


AC n = 


md(x 2 - x- 

= I 


il [V -5 l - r# i *1 + m(z 2 ~ z i ) ^ 2 

[ / q.A F d I J qA F 
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and 


Md(^ - x ; ) r-r- Ji_ +I ,^x S4 ' 

Y I [/ qApd ^ I 


m(y 2 - y^ 2 

qAf 


By iise of equations (l) and (2), 


t = £T - 3 + fa- 

i = % + a + fr 
m' 


If, then, 

°N = °N a + °N,o 


°N " C N,o 
a = — 


% 


a 


. % 
a = -^ L - 


°N, 


a 


C Y - % a P + ^o P 


( C Y ~ C Y,o) 


C N, 


a 


3 - - 


% 


a 


and 


L 


®*Y_P„.p n . p 0d.p *0 

i^d ~ Sx a + °np a 2V P + Si 2V + %° X 


r— Mg 


pd 


) q%d “ Cn tt P + Cn p a 2V a + Cm Q. 2V Cy >° d 
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The last bracketed terms of equations (59) and (40) represent the 
corrections that must be made to ACjj and ACy in order to obtain 

linear plots of Cjj against Cy. Note that the influence of the roll 

can make the correction large even for low roll rates and small ratios 
of I y to I. At high roll rates the correction may become greater than 

the quantity to be determined. Figures 15 and 1^- present experimental 
plots of ACjj aga in st Cjj and ACy against Cy. However, for simplicity 
the points were taken near the oscillation peaks points near 

Cflpeak* points near c Yp eak ^ Where Cjj and Cy are negligible. 

Thus any corrections that must be made to the upper plots of figures 15 
and 14 result entirely from the roll effect on and Z£y. The effect 

of rolling on curves of £Cy against Cy (upper plot, fig. 15) mak es the 
curve appear nonlinear, whereas the effect on ACjj plotted against 
(upper plot, fig. 1^) appears as an increase in scatter. The correction 
for roll given in the preceding equations by use of estimated derivatives 
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reduces the scatter and removes the nonlinearities, as shewn in the lower 
plots of figures 13 and 14. The differences of 10 percent in slopes as 
obtained in the Cjj and Cy planes is not understood. The effect may 

"be real, resulting from different aeroelastic effects on the fins of the 
model, or he within the accuracy of the data. 


The slope s of the curve of 
against Cy is related to cc^, /Vn . 



ACjj plotted against Cjj and ACy 
, \j, and AX by the equation: 



For the test model, the necessary constants are : 


H = 76,600 


2V 

d 


8,120 


J = 2.60 

a 


Iy 

= 0.027 
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Then using the average of slopes from figures 13 and 14 gives 

i 

C Na = 0.0344^125.5 + 0-54 + 0.34) = 4.35 


From equation (13) 


Then 


°^ = 2 (| 2 (lr^ + 0N a.) 

S = «- 5 p^ 2 i( . 0 . 736 ) + 4 . 35 ] 


0^-= 13.5(-13.90 + 4.35) = -129 


Equation (4l) can Ibe rewritten as: 


^ ■ -(I pm 2 ^ + ^ - r^h ** 

i 1 " 21) 


Then 


= -0.0156(125.5 + 0.54 - 0.0023) - 0.00732 


= -I.97- 0.007 = -1.977 
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and 


£Ea _ 5 s& , . hSXL = -0.455 
ac N CH a 4.35 

With Cjj^ and known, the derivative can he obtained from 

equation (^It) and is given by: 



(^ 5 ) 


A value of Cn for the test model could not be determined, since 
Po. 

AA could not be determined with any accuracy at the low roll rates of 
the test. 

It should be emphasized that the biggest drawback in obtaining good 
accuracy is in failure to determine the trim center correctly. In order 
to obtain a reliable trim center, the missile should be disturbed in such 
manner that K is small in magnitude, and the resonance range should be 
avoided. 


CONCLUDING REMARKS 


An analysis has been made of the motions of rolling missiles having 
slight aerodynamic asymmetries. The motions were referred to a body axis 
system. The type of motion encountered by the missile is shown to be 
dependent on the ratio of the rolling velocity to pitching frequency and 
the manner in which the missile is disturbed. The equations of motion were 
found to be 'useful in establishing a te chni que for the reduction of oscil- 
lation data for rolling missiles having body-axis instrumentation. The 
method was applied to experimental results obtained from flight tests of 
a rol li ng missile configuration. 

Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics. 

Langley Field, Va., May 1, 1956. 
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Figure 1.- Examples of rolling trim magnitude for a particular missile configuration. The quan- 


tity K is defined "by equation (2^) . = 0.0392,* 


2 Oat ~ 3} %a 


0 . 
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Figure 2.- Angular orientation of rolling trim for a particular missile. = 0.0392; 

1 t . A rr*_ _ n _ A J _ J 0-1 j A- . _ ± j — - t r\t r \ /-rn .. - . _ J- > - 

=■ ■■ ■ “ = -J; ^Htv, = w. XLH3 tuigj-o jijj in u.cu iiitru uy eQ.uauj.uu uuguoivu 

2(k/ d ) ■ 

values of £a>/oj Q > the sign of is reversed.) 
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pollan, 2 8 4 8 1 

Figure 4.- Typical missile motions below roll resonance for Acd/oIq pos- 
itive. For negative values of the sign of K is reversed. 

Aq = ZA = 0. 





Cycles to 
complete 

pattern,! 4 2 4 I 4 


Figure 5.- Typical missile motions above roll resonance for Au/afc, posi- 
tive. For negative values of Ar/o^>, the sign of K is reversed. 

T^o = => 0. 
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6.- Envelopes of the Cg oscill&tionB for M/Tvq positive 
of £k/’K 0> the sign of K is reversed. 
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For negative values 
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Figure 7.- Response of a rolling symmetrical 
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C n/ N 


^eak alter 
burnout 



Peak before rocket burnout 


■— x Pulse -rocket burning time 


Figure 8.- Amplitude response of a rolling symmetrical missile to the 
disturbance from a pulse rocket having constant thrust. \, = Oj A\ = 

N = \j(d + ~Jj + 4 (d ) 4 (H) f 1 " y) * °R lndica tes resultant 
angle of at tank. 
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